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Abstract. We present high-resolution imaging and low-resolution spectroscopy of the BL Lac object PKS 0537—441 (z = 
0.893) and its environment carried out with the ESO-NTT and VLT telescopes. The observations were designed to clarify, 
whether the properties of PKS 0537—441 are affected by gravitational microlensing due to the claimed detection of a galaxy 
along the line-of-sight to the BL Lac, or whether PKS 0537—441 and its environment act as a lensing system itself, as 
suggested by the detection of several closeby companion galaxies with similar morphologies close to PKS 0537—441. 
Our observations show that neither case seems to be likely. Within our images we did not find a galaxy along the line-of-sight 
to the BL Lac as claimed previously. In addition, our spectroscopy shows that none of the four closest companion galaxies 
(including one new detection by us) is at high redshift. Instead, two of the four nearby companion galaxies to PKS 0537—441 
are within 200 km/s of the systemic velocity of the BL Lac (z = 0.892 and 0.895, respectively). The third companion galaxy 
is at higher redshift (z = 0.947). The fourth companion galaxy shows evidence of Mg II absorption redwards of its systemic 
velocity and is perhaps a mini low-ionization BAL QSO at z = 0.885. If the latter can be confirmed, PKS 0537—441 is the first 
BL Lacertae object being a member of a binary Quasar. 

While we do not find (micro)lensing effects being important for this source or its environment, our observations revealed a 
highly interesting feature. We detected extended [O II] emission in the off-nuclear spectrum of PKS 0537—441, which is most 
likely due to photoionization from the active nucleus, although we can not rule out the possibility that the extended emission is 
due to jet-cloud interaction with the counterjet of PKS 0537—441. 

According to our analysis of the photometric data, PKS 0537—441 seems to be located in a cluster environment as rich as 
Abell type 0-1. This is supported by the detection of four galaxies in the field with similar redshifts as the BL Lac (A z < 
0.002). However, we found serendipitously even more galaxies at somewhat higher redshifts (z = 0.9-1). Thus, PKS 0537—441 
might be located in front of a galaxy cluster at somewhat higher redshift or even be part of a large-scale structure with an 
extension towards the BL Lac. 
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1. Introduction 

BL Lac objects are a relatively rare and extreme subclass of the 
Active Galactic Nuclei (AGN). They are dominated by non- 
thermal emission from radio- up to X-ray frequencies (in some 
cases even up to the 7-ray regime) with mostly weak or even 
absent emission lines in their optical spectra. Within the stan- 
dard AGN picture, their non-thermal radiation emerges from a 
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beamed, relativistic jet closely aligned along the line-of-sight 
to the observer (Urry & Padovani ll995> . 

One of their defining characteristics is their rapid and 
strong variability across the whole electromagnetic spectrum. 
There is common agreement, that the observed variability is in- 
trinsic to the sources (see e.g. the review by Wagner & Witzel 
1995 1, although some BL Lac objects may indeed be micro- 
lensed background QSOs (Ostriker & Vietri ll985t . 

According to the Unified Scheme for AGN, the parent pop- 
ulation of BL Lac objects - i.e. those sources, whose jets are 
misaligned to the observer - are the FR I radio galaxies (Urry & 
Padovani 1995 1. This can be tested by comparing their isotropic 
(unbeamed) properties (e.g. host galaxies or cluster environ- 
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ment), which should be very similar. Host galaxy studies allow 
us to test the micro-lensing scenario, since here an offset be- 
tween the center of a detected galaxy and the AGN itself could 
be present. In this case the galaxy is presumably a system along 
the line of sight rather than the host of the AGN. 

In recent years, in more than 100 BL Lac objects with z < 
0.5 a host galaxy has been detected (e.g. Falomo 1996 , Wurtz et 
al. [T996l Heidt et al. fl999l Falomo & Kotilainen [T999l Scarpa 
et al. 120001 Nilsson et al. 2003 ) and there is common agree- 
ment that BL Lac host galaxies are luminous, giant elliptical 
galaxies. Host galaxy data for BL Lacs at higher redshifts are 
scarce, so any evolutionary trends can not be investigated. The 
cluster environment has been studied by Wurtz et al. ( 1997 ) for 
a sample of ~ 50 sources up to z = 0.65 (mostly at z < 0.4). 
They found an apparent evolution towards denser cluster en- 
vironments with redshift. However, the apparent evolution is 
tied to observations of a few BL Lac objects between z = 0.5 
and 0.65 in their sample only. Similar conclusions have been 
drawn by Fried et al. ( 1993 1 from a study of the 1 Jy sample of 
BL Lacs (Stickel et al. 1993 ). Unfortunately, due to the limited 
depth of their observations, they were not able to evaluate the 
cluster environment for high-redshift (z = 0.8-1) 1 Jy BL Lacs 
(including PKS 0537-441). 

Up to now, there is only one clear example of a lensed BL 
Lac (B2 0218+35.7, Patnaik et al. FT993l . Other objects, whose 
properties may be affected by gravitational (micro)lensing are 
AO 235+164, PKS 0537-441 and B2 1308+326 (see Heidt 
I1999l for details). These three BL Lacs all have similar proper- 
ties. They are at relatively high redshift (z ~ 1), are among the 
most luminous BL Lac objects (Mr < -27.5) and have histor- 
ically shown the largest variability amplitudes (A m > 5 mag, 
see Tab. 1 in Stickel et al. [1993 1. A few more candidates de- 
tected during the course of the BL Lac HST snap survey are 
discussed in Scarpa et al. (119991 . 

Of particular interest is the BL Lac object PKS 0537-441 
(z = 0.892, Lewis & Ibata 2000 1, which is a member of the 1 
Jy sample of BL Lac objects (Stickel et al. 1993 ). It displayed 
intraday variability at radio (Romero et al. 119951 and optical 
frequencies (Heidt & Wagner 1996 Romero et al.|2000l, and 
is one of the few BL Lac objects, where 7-ray emission has 
been detected (Hartman et al. 1999 1. Its continuum variations 
throughout the entire electromagnetic spectrum have been dis- 
cussed by Pian et al. ( 2002 1. Interestingly, Romero et al. ( 1995 1 
derived from their radio variability measurements brightness 
temperatures in excess of 10 21 K, which is 9 orders of magni- 
tude above the inverse Compton limit. Based on their claimed 
detection of a galaxy along the line-of-sight to PKS 0537—441 
Stickel et al. (1988 see below) discussed, whether the BL 
Lac could be microlensed by small masses in this foreground 
galaxy. On the other hand, spectral-index variations found from 
optical two-color photometry by Romero et al. ( 2000) seem to 
rule out the microlensing hypothesis at least for these observa- 
tions. 

PKS 0537—441 was also subject to several broad-band 
imaging studies through a R-filter, however, with puzzling re- 
sults. Stickel et al. (1988 1 claimed the detection of a fore- 
ground disk galaxy to PKS 0537—441, with a morphology 
and brightness similar to a galaxy 11 "east of the BL Lac 



at z = 0.186. They concluded that the observed properties 
of PKS 0537—441 are influenced by gravitational lensing by 
the foreground galaxy. In subsequent studies Falomo et al. 
(I1992t and Kotilainen et al. dl998> did not detect the fore- 
ground galaxy reported by Stickel et al. (1988) and found 
PKS 0537—441 completely unresolved thus weakening the 
lensing hypothesis. The former noted two faint companion 
galaxies close to the BL Lac at ~ 4"distance. The BL Lac ob- 
ject also appeared unresolved in HST observations (Scarpa et 
al. |2)00). In a reanalysis of the same HST-image mentioned 
above, Lewis & Ibata (2000) found indications for a galac- 
tic component after the subtraction of a scaled PSF (which 
was recently shown to be due to scattered light by Pian et al. 
2002). Additionally, they detected 3 companion galaxies close 
to PKS 0537-441 (2 of which had already been found by 
Falomo et al. 1992) and discussed an alternative scenario in 
which PKS 0537—441 is not subject to lensing effects, but is 
rather a member of a group of galaxies which lense a distant 
background galaxy. As a consequence at least two of the three 
companion galaxies could be magnified and split images of a 
distant background source at high redshift. 

Obviously, the results and the interpretation of the results 
for PKS 0537—441 differ grossly. One major drawback of these 
studies is that the images were either taken under moderate see- 
ing conditions or suffer from relatively short exposure times. 
Except for PKS 0537-441 itself and the nearby galaxy ll"to 
the east, no spectroscopy in the field of PKS 0537—441 has 
been carried out until now. In order to clarify the situation, we 
have undertaken a very deep imaging and multi-object spectro- 
scopic study of PKS 0537—441 and its environment using the 
NTT and VLT telescopes at ESO, the results of which are pre- 
sented here. While we will show that lensing effects are most 
likely not important in this source, several highly interesting 
features are observed. Amongst them are the presence of ex- 
tended [O II] emission on the opposite side of the radio jet in 
PKS 0537-441 and the possible detection of a binary QSO. 

To facilitate comparison with other studies, Ho = 50 km 
s _1 Mpc -1 and qo = is assumed. 

2. Observations and data reduction 

2. 1. Imaging 

A very deep I-band image of PKS 0537—441 was acquired at 
the NTT in service mode on the night November 7/8 1999. The 
use of the I-filter instead of a R-filter (which was used in pre- 
vious studies) has the advantage, that it samples the spectrum 
redwards of the 4000A break at z = 0.893 thus increasing the 
chance to detect the host gfalaxy of PKS 0537—441. We used 
SUSI2, which consists of a mosaic of two 2kx4k EEV-CCDs. 
The frames were binned by a factor of two, which gives a scale 
of C.'162/pixel. Thus the field of view is ~ 5. '5 x 5 '.5. During 
the photometric night 43 images (3x1 min, 40 x 3 min) 
resulting in an integration time of 123 min in total were taken. 
PKS 0537—441 was placed close to the center of one of the two 
CCDs (CCD #45, which covers the eastern part of the mosaic). 
The short individual images were necessary to avoid saturation 
of the BL Lac nucleus. In order to use the frames for fringe 
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correction and superflat purposes, jittered images were taken. 
We used a random walk jitter pattern within a rectangular box 
of 30" border length centered on the central position. During 
the night standard stars from Landolt ( 1992) were observed to 
set the zero point. 

The individual images were first corrected for bias and then 
for pixel-to-pixel variations. For the latter we used twilight flat- 
fields. A superflat was created from the 43 bias-subtracted and 
flatfielded images to correct for large-scale gradients. To re- 
move the fringes a scaled fringe template derived from I-band 
images taken during this night was subtracted. Finally, the im- 
ages were cleaned of cosmic ray hits, aligned and summed. The 
resulting image has a FWHM of (X'85. 

Additionally, a short R-band image of PKS 0537-441 has 
been taken with the VLT-UT1 and FORS on the night October, 
10-11 1999 for the preparation of the spectroscopic observa- 
tions. The data reduction was similar as for the I-band images. 



2.2. Spectroscopy 

Low-resolution spectra of ~ 46 objects in the field of 
PKS 0537-441 were collected with FORS1 at the VLT1 
(Antu) during two clear nights on November 19-21, 2000. 
Observations were done in MOS (Multi-Object Spectroscopy) 
mode, in which FORS1 provides 19 slitlets of ~ 22"length 
in a vertical array. Two MOS-setups were designed. Setup A 
was observed for 1 hour (3 x 20 min) at PA = 49? 6 (counted 
counterclockwise from N), whereas setup B was observed for 
2 hours (3 x 40 min) at PA = 32?6. We used the grism 
1501, which gave us a spectral scale of ~ 5.5 A/pixel. The 
slit width was set to 1". The wavelength range covered was 
~ 4000 - 10000A. Slitlet #9 for setup A covered two brighter 
nearby companion objects of PKS 0537—441 (labeled Al and 
A2 in Fig. whereas the slitlet #9 for setup B covered 
two fainter nearby companion galaxies of PKS 0537—441 (la- 
beled Bl and B2 in Fig.0. The remaining slitlets were filled 
serendipitously with galaxies in the field in order to search for 
galaxies at the redshift of PKS 0537—441. In some cases we 
were able to put two objects onto one slitlet. The mean DIMM 
seeing was 0?.'6 FWHM for setup A and 1'.'05 FWHM for setup 
B. Finally, a short (lOmin) long-slit spectrum across the nu- 
cleus of PKS 0537-441 at PA = with the same slit width 
and grism was observed on the night November, 21-22 during 
a period of poor seeing conditions (1'.'8 FWHM). At the end 
of each of the nights, spectrophotometric standards from Oke 
J1990i were observed. 

The data reduction of the individual spectra (bias subtrac- 
tion, flatfielding, cosmic ray removal, sky subtraction, wave- 
length calibration, etc.) was performed using standard MIDAS 
routines. The individual spectra of each slitlet were averaged 
and smoothed using a running mean with a width of 3 pixels. 
The FWHM spectral resolution measured from strong night sky 
emission lines is ~ 25A. 



3. Analysis and results 

3. 1. Imaging 

3.1 .1 . The host galaxy and nearby environment 

In order to search for the host galaxy or a galaxy along the line- 
of-sight to PKS 0537—441 we applied a fully 2— dimensional 
fitting procedure to the image (see Nilsson et al. 1999 and Heidt 
et al. 1 19991 for details). We first extracted a PSF by averaging 
the three brightest stars in the field of PKS 0537-441. Only 
stars present on the same CCD as PKS 0537—441 were used. 
Unfortunately, PKS 0537—441 was bright during the observa- 
tions. Thus, even the brightest stars in the field usable for the 
extraction of the PSF were at least 1.5 mag fainter than the BL 
Lac itself. In order to improve the S/N in the outer parts of 
the PSF, we smoothed the PSF outside the radius of 2'.'4 (~ 3 
FWHM) using a Gaussian kernel. In the inner part (r < 2'.'4) 
no smoothing was applied in order to accurately preserve the 
rapidly changing brightness profile close to the center of the 
PSF. This composite PSF was then used for the model fits. 

We fitted three different models to the BL Lac. One rep- 
resenting the AGN (scaled PSF), one representing an AGN 
+ bulge (scaled PSF + convolved de Vaucouleurs model) and 
one representing an AGN + disk (scaled PSF + convolved disk 
model). We did not allow for any offset between the AGN and 
model galaxy, the position angle and ellipticity of the putative 
host galaxy were set to zero. Prior to fitting, we masked care- 
fully all regions affected by nearby companion objects. 

Already the fit with a scaled PSF was an adequate descrip- 
tion of the observed light distribution of PKS 0537-441, i.e. 
we could neither detect the host galaxy of PKS 0537—441 
nor a potential galaxy along the line of sight to the BL Lac. 
PKS 0537—441 appears completely unresolved. Its brightness 
was mi = 14.45, which is relatively bright for this source. The 
image of PKS 0537—441 and its residuals after subtraction 
of the best fit AGN is shown in Fig.^ the 1-dim profiles of 
PKS 0537-441 and a scaled PSF in Fig.El 

As can be seen from Fig. [fl there are residuals in 
the center of PKS 0537-441, which are typical for PSF- 
subtraction. Remarkably, our PSF traces the surface brightness 
of PKS 0537-441 out to a radius of about 8" covering -12 
magnitudes in brightness (Fig. |2ji. There are some slight dif- 
ferences at about 2-3" from the core (seen also as remaining 
light to the north and west of the BL Lac in Fig. 0, whereas 
at radii > 8" the PSF clearly underestimates the flux from 
PKS 0537-441. Both features are due to to the imperfect PSF. 

We have estimated upper limits for the host galaxy bright- 
ness using simulated images of PKS 0537—441. The simulated 
images consisted of a core component with mi = 14.45 and 
a host galaxy with effective radius r c ff = 0"95 (10 kpc at z 
= 0.893). Using progressively fainter host galaxies we deter- 
mined the highest host magnitude that allowed us to make a 
host detection. For the host galaxy to be detected we required 
r"g > 5cr rc ff , where r°g is the average effective radius and cr re ff 
is the rms scatter in a set of 30 simulations. Using this method 
we estimate that a host galaxy with mi < 19.0 would have been 
detected. 
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Fig. 1. Left) NTT+SUSI2 I-band image of PKS 0537-441. FOV is 30" x 30". Nort h is up, east to the left, FWHM ~ Of! 85. The 
galaxy Gl 11" to the east of PKS 0537-441 at z = 0.186 is discussed in Stickel et al. J1988l >. Right) Same image after subtraction 
of a scaled PSF. The four closest companion galaxies to PKS 0537—441 are labeled. Labeling for Al, A2 and Bl follows Falomo 
et al. ( 19921 and Lewis & Ibata (2000 1. B2 is a new detection by us. The orientation of the MOS setups A and B as well as of the 
longslit spectrum is also indicated. 



PKS 0537-441 




radius/[arcsec] 

Fig. 2. Surface brightness profile of PKS 0537-441 (asterisks) 
compared to a scaled PSF (solid line). 



Inspection of Fig.^reveals four resolved galaxies within 
7"of PKS 0537-441. The three brightest ones ( Al, A2 and 
Bl) have already been detected by Falomo et al. (0992 1 and 
Lewis & Ibata (2000 1. In addition, there is a fourth companion 
galaxy - 6'.'5 to the north of PKS 0537-441 (labeled B2 in 
Fig. [D, which is also barely visible in Fig. 1 of Falomo et al. 
dl992i . The brightness of these objects is mi ^21.0, 21.2, 22.0 
and 23.0 and their projected distances from the BL Lac are 3"9, 
3"8, 5"and 6"4 for Al, A2, Bl and B2, respectively. 



3.1 .2. Cluster environment 

In order to analyze the cluster environment of PKS 0537—441 
we first created an object catalog with positions, magni- 
tudes (MAG_BEST) and a star/galaxy classification index 
(CLASS.STAR) using SExtractor (Bertin & Arnouts [T996V 
MAG_BEST is an isophotal magnitude which considers de- 
blending of objects and applies an aperture flux correction at 
the edges of sources. The CLASS.STAR parameter can be in 
the range between (galaxy) and 1 (point source). This value 
is calculated by a comparison of the FWHM of an image and 
the morphology of objects. It considers both the ellipticity of 
objects as well as their ratio of maximum brightness to their 
extension. 

Within the field of 2 (CCDs) x 2. '32 x 5 .'15 we detected 
~ 1200 sources. After a careful visual inspection of the detec- 
tions, ~ 60 obviously wrong entries were removed (false de- 
tections close to the edges of the frame or imperfectly removed 
cosmic ray hits) and a few sources were added (which were not 
detected by SExtractor but clearly present). In the latter case we 
performed aperture photometry with standard Midas routines. 
The completeness limit derived from the peak in the differen- 
tial galaxy number counts is mi ~ 24.5. Finally, the projected 
distance (in kpc) of each source to the BL Lac was calculated. 

To determine the strength of galaxy clustering around 
PKS 0537—441 we measured the surface galaxy density as a 
function of projected distance from the BL Lac. All sources 
within five circular annuli with equidistant radii from 100 kpc 
to 500 kpc centered on PKS 0537—441 were counted and their 
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Radius 0(in kpc) 

Fig. 3. Left) The field around PKS 0537—441 from the NTT I-band image. The 500 kpc environment is indicated by a circle 
centered on PKS 0537—441. Plotted are all objects with CLASS_STAR<0.5 (i.e. galaxies). Galaxies in the magnitude range 
m g < m < m g + 3 are marked by a square. PKS 0537—441 is included as a galaxy, according to the magnitude-redshift relation 
by Eales ( 1985) we used m g (I) = 19.5. A density excess towards the BL Lac can easily be seen. There are 24 galaxies within and 
47 galaxies outside the 500kpc region, i.e. 1/3 of the galaxies are within 1/6 of the area covered by this CCD. Right) The galaxy 
density excess u>{ff) as a function of radial distance to the BL Lac. The mean field galaxy density is indicated by the horizontal 
line. It has been derived from the entire field (i.e. both 2kx4k CCDs) excluding the 500kpc environment of the BL Lac. The top 
panel shows u>{6) for all galaxies, whereas in the bottom panel only galaxies within the range m g < m < m g + 3 have been 
used. For both diagrams only galaxies (CLASS_STAR<0.5) are considered. 



resulting counts compared to those derived for field sources. 
Precisely, we calculated 



w(0 + AO) 



N B ( 



Ad) 



N F (9 + A9) 
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where Nb (9 + A9) are the number counts within a given area 
9 + A6 from the BL Lac and N F {9 + A6) are the counts which 
would be expected if the galaxies are distributed as in the field. 
The mean field galaxy density N / arcmin 2 was calculated 
from galaxy counts at 9 > 500 kpc from which N F (9 + A6) 
has been scaled according to the area covered. Both, cumula- 
tive counts as well as differential counts have been measured. 

Fig- E (top right) shows uj{9 + A9) by using all galaxies 
detected (CLASSJSTAR < 0.5). Even without any restriction 
of the galaxies to a certain magnitude range a density excess 
towards PKS 0537—441 can be seen. The clustering signal be- 
comes even stronger if we use only galaxies in a magnitude 
range, which has the highest probability to be at the BL Lac 
redshift (Fig. [5] bottom right). For the definition of the magni- 



tude range we followed Hill and Lilly (1991 hereafter HL91) 
and considered only galaxies within m g < m < m g + 3. The 
typical apparent magnitude m g of a giant elliptical galaxy at a 
given redshift was obtained from the magnitude-redshift rela- 



tion for radio galaxies of Eales ( 1985 1 
m g (R) = 21.05 + 5.3logz. 

Since we observed in the I-band, m g (I) was estimated by sub- 
tracting the R — I color of an early type galaxy (k-correction 
and galactic extinction towards PKS 0537—441 had also been 
taken into account). 

For an estimate of the cluster richness and a direct com- 
parison to other studies, we calculated the excess of galaxy 
counts N 5 as defined by HL91. This is the number of galax- 
ies ranging from m g < m < m g + 3 within the 0.5 Mpc 
radius around the BL Lac minus the average number counts 
of the field galaxy distribution. For PKS 0537—441 we find 
iVo.5 = 11.5 ±3.5. This corresponds to a spatial cross- 
correlation amplitude of the order of B qg « 435 following the 
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relation B qg = (37.8 ± 10.9)7V .5 given by Wold et al. <2000l, 
and hence to a cluster of Abell richness class 0-1 (if compared 
to the calibration in Tab. 4 and Fig. 3 of HL91). 

3.2. Spectroscopy 

Preliminary results of the spectroscopic observations have been 
described by us in Heidt et al. OOOlal 1200 lbl 120031 . These 
were solely based on the identification of emission/absorption 
lines in the spectra. For the present paper more accurate red- 
shifts could be determined by a more sophisticated analysis, 
which further allowed us to roughly classify the objects into 
galaxy types. 

The spectra were first corrected for the atmospheric A-band 
and B-band using the standard stars during the same nights. 
This is very important, since the A-band coincides at z ~ 0.9 
(the BL Lacs redshift) with the 4000 A break in early-type 
galaxies. Then we searched for obvious emission/absorption 
line features in our spectra and determined their centroids by 
fitting Gaussians to them. In order to avoid misidentifications 
due to artefacts (e.g. imperfect removed cosmic ray signatures), 
all the 20 min and 40 min spectra were inspected and com- 
pared individually. Finally, we fitted interactively five different 
spectral energy distributions (SEDs) to the spectra. As start- 
ing values we used preliminary redshifts obtained from emis- 
sion/absorption lines. The five different SEDs used for the fit- 
ting procedure were derived from spectroscopic observations of 
about 500 galaxies in the FORS Deep Field 1 (S. Noll, in prep.). 
They have the advantage that they were derived with the same 
spectroscopic setup, instrument and telescope as the observa- 
tions presented here. The galaxy types covered are: SED1 = 
E-S0, SED2 = Sa-Sc, SED3 = Sd/Irr, SED4 = starburst, SED5 
= extreme starburst. Four of the SED5 sources (#a6, #a7, #b6 
and #bl0) may be AGN. They all have the commonly used flux 
ratios [O III] 5007/[O II] > 1, [O III] 4959/H/? > 1 and [O III] 
5007/H/3 > 3 for a classification as an AGN. Unfortunately, 
due to the relatively high redshift of the sources, the impor- 
tant diagnostic lines [N II] 6583/Ha are either in the extreme 
red part of the spectrum dominated by night sky emission or 
outside the spectral range covered. The ratio [O III] 5007/H/3 
vs. [ N II] 6583/Ha is a robust diagnostic diagram to separate 
starburst galaxies from AGN (Veilleux & Osterbrock l 1987b . 

The approach of fitting SEDs is especially helpful for the 
redshift determination for galaxies dominated by an old stel- 
lar population. The error in the redshift determination is Az 
~ 0.001 for sources with at least 2 emission lines and Az < 
0.003 for the remaining sources. The latter error was derived 
by slightly varying the redshift of the determined SEDs for in- 
dividual galaxies until an obvious mismatch was apparent. 

The results are summarized in Tab. ^ Column 1 identifies 
the slitlet, where the objects were located, columns 2 and 3 
give the offsets with respect to PKS 0537—441 followed by 
the apparent and absolute R and I magnitudes in columns 4 
to 7. The apparent R and I magnitudes were measured using 
SExtractor (Berlin & ArnoutsQl96j. Since PKS 0537-441 is 
strongly saturated on the R-band acquisition image taken for 




see Heidt et al. 1 2003 l for details on the FORS Deep Field 



Fig. 5. Distribution of the objects from the 2 MOS-setups 
with redshift and type. Hatched = SED1-3, double hatched = 
SED4/5 and black = AGN, respectively. 



the preparation of the MOS setups, we do not give R magni- 
tudes for the companion galaxies Al, A2, Bl and B2. The R- 
and I-images do not fully overlap. Hence, I magnitudes could 
not be derived for all sources. No k-correction was applied to 
the absolute magnitudes, because the classification into SEDs 
is very rough. Column 8 gives the redshift and column 9 the 
rough SED classification. Finally, in column 10, the detected 
lines are listed and comments are given. In Fig. |4] representa- 
tive spectra of objects with SED 1-5 and one AGN candidate 
are shown. 

For 36 out of the 46 spectra in total, we were able to deter- 
mine the nature of the sources. Their R and I magnitudes range 
from 18 to 24.5, most have hir > 22.5. Their absolute magni- 
tudes range from -15.3 to -22.6 (median -20.6) in R and from - 
16 to -30 (median -21.9) in I. Remarkable is the extreme bright- 
ness of the BL Lac (-30.0) and brightness of the two nearby 
companions Al (-23.5, QSO) and A2 (-23.5, early-type at z = 
0.947). We found 32 galaxies (4 SED1, 5 SED2, 12 SED3, 4 
SED4 and 7 SED5), 2 AGN (PKS 0537-441 and Al, which 
will be discussed in the next section) and 2 stars. 

For about 60% (19/34) of the galaxies, at least two emis- 
sion/absorption lines could be detected, which allowed a se- 
cure redshift measurement. In 14 sources we could only detect 
a single narrow emission line, which we identified with [O II] 
in all cases. This identification was supported in most cases by 
either a jump of the continuum redwards of the [O II] line indi- 
cating the 4000 A break or a blue continuum typical for strong 
star forming galaxies. Finally, in companion galaxy Al, we de- 
tected a broad emission line which we identified as Mg II. 

The redshifts of the objects range from 0.07 up to 1.3, 19 
are at z > 0.8. Most objects are dominated by [O II] emission, 
implying relatively high star formation rates. We found also 
several early-type galaxies at z > 0.9. Four objects have red- 
shifts similar to PKS 0537-441 (A z < 0.002). The distribu- 
tion of objects with redshift and type is shown in Fig.|5] 

3.2.1 . PKS 0537-441 and A1 , A2, B1 & B2 

In Fig. |6] the 1-dim spectra of PKS 0537-441 and the com- 
panion galaxies Al, A2, Bl and B2 are displayed. Although 
the companion galaxies are about 6-8 mag fainter than the 
BL Lac, their spectra were sufficiently well separated from 
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Those objects were labeled *a and *b, respectively. 



PKS 0537—441 to measure not only characteristic lines, but 
also to detect the continuum. The spectra are discussed below. 

PKS 0537-441: In the off-nuclear spectrum of 
PKS 0537-441, Mg II at 5300 A and [O II] at 7052 A 
are detected, which give z = 0.893 on average. In the spectrum 
across the nucleus, only Mg II at 5300 A (z = 0.894), but not 
[O II] was unambiguously detected. Both lines have previously 
been found by Lewis & Ibata (2000 1 on their nuclear spectra 
with higher resolution. They give z = 0.892 in good agreement 
with our measurement. We also note that [O II] has not been 



detected by Stickel et al. (119931 and by Falomo et al. ( 1994 1, 
although it may be present in the spectrum shown by Stickel et 
al. (cf Fig. 2 in their paper). 

Al: The blue spectrum of this object shows a very broad 
asymmetric line at 5275 A and a narrow line at 7052 A, which 
are most likely Mg II and [O II]. The Mg II line is ~ 25 A 
blueshifted with respect to the Mg II line in PKS 0537-441 
and gives z = 0.885. It can not be Lya at z = 3.338 since the 
continuum extends towards the blue and no obvious Ly-break 
is present. We emphasize here that the broad asymmetric line 
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Fig. 4. Spectra of representative galaxies and of an AGN candidate detected in the field of PKS 0537— 441 . Templates of average 
spectral energy distributions of galaxies of different Hubble type (SED1-5) are overplotted on the observed galaxies' spectra. 
The error spectrum (dashed line) is displayed as well. 



at 5275 A can also not be scattered light from the Mg II line in 
PKS 0537-441 due to reflections inside FORS. This line in Al 
can always be seen with a similar structure in all 3 individual 
images. Moreover, reflections inside FORS normally scale as 
10 -4 (W. Seifert, priv. com.), whereas we observe here a dif- 
ference of the Mg II lines of the order of 10 _1 . However, the 
center of the [O II] line in Al is at the same position as the [O 
II] line in the off-nuclear spectrum of PKS 0537—441. 

In order to investigate the apparent discrepancy between 
the Mg II and [O II] lines, we subtracted the continuum on the 
2-dim spectrum of Al, A2 and PKS 0537-441. This was ob- 
tained by fitting and subtracting row-by-row a smoothed cu- 
bic spline from the 2-dim spectrum. Great care was taken to 
match the continuum next to 5300 A and 7050 A as well as 
possible. In Fig.[7]the 2-dim averaged, smoothed spectrum and 
the resulting continuum-subtracted spectrum are shown. There 
are several remarkable features. First of all, the Mg II line of 
PKS 0537-441 and the Mg II line of Al are well separated. It 
is also obvious, that the Mg II line of Al is very broad and its 
center somewhat blueshifted with respect to the Mg II line in 
the BL Lac. By contrast, the [O II] line of PKS 0537-441 ex- 
tends from the central position (at about 1'.'7 from the core) 
with decreasing intensity without any blueshift towards Al. 



The projected length of that feature is about 4"(~ 40 kpc at 
z = 0.893). There might also be some (kinematic) substruc- 
ture (wiggles) present, but this could be due to the continuum- 
subtraction procedure. Since the central position of Mg II and 
[O II] in the BL Lac give the same redshift within the errors 
contrary to the central position of Mg II in A 1 and since the 
intensity of the [O II] line decreases towards Al without any 
blueshift, we believe that [O II] is intrinsic to PKS 0537-441. 
Curiously, [O II] does not extend in the opposite direction. A 
physically plausible explanation for this complex configuration 
will be discussed in the next section. 

A2: The spectrum of this galaxy is dominated by a red con- 
tinuum showing Ca K+H and Mg b absorption lines character- 
istic for early-type galaxies with SED2. No emission lines are 
detected. The redshift of this galaxy is z = 0.947. 

B 1 : This galaxy has also a relatively red spectrum (SED2) 
showing Ca K+H and the g-band in absorption as well as a faint 
[O II] emission line. The resulting redshift is z = 0.895, close 
to the redshift of PKS 0537-441 . 

B2: Similarly to Al, the spectrum of this source has a 
blue continuum (SED5) and a relatively strong narrow emis- 
sion line, which we identify with [O II]. The resulting redshift 
is z = 0.892, again, close to the redshift of PKS 0537-441. 
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Fig. 6. Spectra of PKS 0537—441 (off-nucleus from setup A and across the core from longslit observations) and of the 4 com- 
panion galaxies Al, A2, Bl and B2. For the companion galaxies the error spectrum (dashed line) is included. 



4. Discussion 

Although our I-band image of PKS 0537—441 is probably the 
deepest ever taken of that field, we were not able to resolve the 
host galaxy of the BL Lac. This might be at least in part due to 
the absence of sufficiently bright stars in the field which pre- 
vented us to extract a very good PSF. The lower limit to the 
host galaxy magnitude based on our simulations is mi = 19.0, 
which would result in Mi = -26.8 (k-correction included). With 
a typical color of ellipticals of R-I = 0.7 (Fukugita et al. 1995 1 
this would give Mr = -26.1. Thus, using a typical BL Lac host 
magnitude of Mr = -23.5, our observations should have been at 
least 3 mag deeper given the less than optimal PSF. On the other 
hand, we can rule out an intervening galaxy as bright as the 
early-type companion galaxy Gl (mi ~ 17.4) 11" to the east, 
which was seen by Stickel et al. ( 1988 1. This confirms the ob- 
servations by Falomo et al. ( 1992 1, Kotilainen et al. ( 1998 1 and 
Scarpa et al. (2000 1. Clearly, we can not rule out an even fainter 
galaxy along the line of sight (our detection limit is about L* 
at z ~ 0.2). Nevertheless, there is no strong evidence that the 
properties of this very luminous and strongly variable BL Lac 
nucleus are affected by gravitational microlensing effects. This 
is further supported by the absence of absorption lines in HST 
FOS spectra (Bechtold et al. 2002 1. A deep high-resolution im- 



age of PKS 0537—441 in a very low state is required to unam- 
biguously detect the host galaxy or a faint galaxy along the line 
of sight, but this is hard to accomplish. 

To test, if some of the nearby companion galaxies are 
lensed background sources, we have carried out multi-object 
spectroscopy in the field of PKS 0537—441 including the four 
companion objects. None of the four nearby companion galax- 
ies was found to be at high redshift, i.e. it also does not seem 
that the environment of PKS 0537—441 acts as a lensing sys- 
tem. 

On the other hand, our spectra show that PKS 0537-441 
and the companion galaxies Al, A2, Bl & B2 form a very com- 
plex system. Whereas A2 is an early-type galaxy at z = 0.947, 
the companion galaxies Bl and B2 (z = 0.895 and 0.892, re- 
spectively) are within 200 km/s systemic velocity of the BL 
Lac, i.e. they are very close to PKS 0537—441. Due to their 
proximity to PKS 0537—441 they could trigger the strong ac- 
tivity in this source via gravitational interaction. 

The relation between PKS 0537—441 and Al is very hard 
to interpret. Whereas we detected Mg II in both, the nuclear 
and off-nuclear spectrum of PKS 0537—441, [O II] was only 
detected in the latter. If one compares the relative strengths of 
Mg II in the two spectra (cf Fig. |6j, one would also expect 
to see [O II] in the nuclear spectrum, which is, however, not 
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Fig. 7. Top) 2-dim spectrum of slitlet #9 from setup A. The positions of the companion galaxies Al and A2 are displayed. A 
cosmic ray hit (C.) is also indicated. The spectra of both companion galaxies are well separated from the off-nuclear component 
of PKS 0537— 441 . Note the relatively blue continuum of Al and the relatively red continuum of A2. The broad Mg II emission at 
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spectrum after subtraction of the continuum. Note that the Mg II emission from Al is well separated from the Mg II emission of 
PKS 0537—441, whereas the [O II] emission extends from the BL Lac towards Al. Note also that the [O II] emission does not 
extend towards the opposite direction. See text for details. 



present. We also failed to detect [O II] in the nuclear 2-dim 
spectrum by subtracting the continuum in the same way as de- 
scribed for the off-nuclear spectrum. This points to the fact that 
the [O II] emission is only excited in that region, which is sup- 
ported by the detection of extended, asymmetric [O II] emission 
with decreasing intensity towards Al in the 2-dim off-nuclear 
spectrum. 

There are several mechanisms, which may give rise to the 
[O II] emission starting ~ 20 kpc from the core and extend- 
ing ~ 40 kpc outwards: locally induced star formation due 
to gravitational interaction, photoionization by the central ac- 
tive nucleus or shock-ionization via interaction between radio- 
emitting components and the ambient gas (jet-cloud interac- 
tion). 

In order to induce star formation by gravitational inter- 
action a companion galaxy or at least a merger remnant is 
required. We think that this is unlikely here. No companion 
galaxy is detected in that region (Al is at lower redshift) and 
no obvious tidal features are present at the location of the ex- 
tended [O II] region. Moreover, elliptical galaxies normally do 
not possess much gas at ~ 3 r (under the assumption that the 
host galaxy of PKS 0537—441 is a typical host galaxy with r c 
= 10 kpc). 

Whether photoionization by the active nucleus or jet-cloud 
interaction is the dominant mechanism depends mainly on the 
size and extent of the radio sources with respect to the emission 
line regions (see review by Tadhunter 2002 1. Photoionization 
is normally the dominant mechanism for those AGN, where 



the radio source extends well beyond the emission line regions 
(Tadhunter 2002 1 although kinematical evidence for jet-cloud 
interaction has recently been found in exactly such sources (e.g. 
Solorzano-Inarrea et al. 2001 1. 

With our limited information available, we are not able 
to distinguish between or set even constraints on the last two 
possibilities described above. Whereas photoionization by the 
active nucleus is certainly the more realistic explanation (ex- 
tended line-emission up to 20 kpc from the nucleus has been 
detected e.g. in the BL Lac PKS 0521-36, Boisson et al. H989l . 
the alternative can not be ruled out. PKS 0537—441 is one of 
the few BL Lacs classified as a FR II radio source (Rector & 
Stocke 2001 and references therein) and shows in radio images 
a curved jet-like structure leading to the west (Cassaro et al. 
119991 . Remarkably, our extended [O II] emission is present al- 
most exactly on the opposite side of the jet of PKS 0537—441 
extending ~ 40 kpc in length. This is illustrated in Fig. [8] 
where we show a comparison of the location of the jet in 
PKS 0537—441 and the extended [O II] emission. We may see 
here at least some signature of a counterjet in this BL Lac. 
Evidence for a counterjet has been detected via kinematic sig- 
natures of jet-cloud interaction (line-splitting in [O III]) in 3C 
120 (Axon et al. fT989l . We note, that some kinematic signa- 
tures in the extended [O II] emission in our 2-dim off-nuclear 
spectrum of PKS 0537—441 may be present, although this may 
be due to the continuum-subtraction procedure. 

Further inspection of Fig. [8] reveals that the jet in 
PKS 0537—441 apparently turns over close to BL This may 
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Fig. 8. Same as Fig. □ with the VLA 1360 MHz image of 
PKS 0537—441 superimposed (courtesy of P. Cassaro). The 
region, where the extended off-nuclear [O II] emission towards 
Al has been detected is indicated. It is obviously on the oppo- 
site side to the jet. Note also the coincidence of the location of 
Bl and the turnover of the jet. FOV is 30" x 30", north is up, 
east to the left. 

be due to interaction of the jet with Bl, where a pressure gradi- 
ent deflects the jet considerably. Again, whether this is correct 
or not can not be tested by our observations but it is not un- 
usual. A similar feature has been observed in the low-redshift 
BL Lac 3C 371 albeit on smaller spatial scales (Wrobel & Lind 
[T990I Nilsson et al. [T997l . 

The nature of Al is also unclear. As we have argued in 
the previous section, the emission line with some absorption 
towards the red is most likely Mg II at z = 0.885. With M T 
= -23.5 we classify Al therefore as a QSO (note that our I- 
filter samples almost exactly the B-band at z = 0.885). Fitting a 
Gaussian to the absorbed emission line gives a width of the line 
of - 100A FWHM (~ 5000 km/s). Corrected for redshift this 
gives a FWHM of ~ 53A, which is is much larger as typically 
measured in QSOs (~ 35A FWHM, Vanden Berk et al. l200Tl 
If the absorption is due to an absorbing system along the line of 
sight to Al one would also expect absorption in the spectrum 
of PKS 0537—441, which is not observed. Therefore, the ab- 
sorption must be intrinsic to the source and Al could hence be 
a mini low-ionization BAL QSO (LoBAL). 

Although classical BAL QSOs have their absorption blue- 
wards of their systemic redshifts, a few BAL QSOs with ab- 
sorption redwards have been detected in the SDSS (Hall et al. 
2002 1. The absorption redwards of the systemic redshift can 
be either due to infalling gas crossing our line of sight or due 
to a rotation-dominated disk wind (Hall et al. 2002 1. Since we 
do not detect unambiguously [O II] at ~ 7025 A (z = 0.885) 
in Al indicative of a relatively high covering factor, infall of 
gas could be the dominating absorption mechanism. There is 
at least indirect support for the identification of A 1 as a mini 
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Fig. 9. Comparison of the Mg II lines in Al (solid line) and 
the BL Lac (long dashed line). The latter was shifted by 25 A 
towards the blue and scaled up by a factor of 2. The contin- 
uum was subtracted in both spectra. Note the good match of 
the wings of the lines. For comparison the original position of 
the center of the Mg II line in PKS 0537-441 is indicated. 



LoBAL by the non-detection of [O II], which is often absent 
or weak in Mg II LoBALs (Boroson & Meyers 1992 1 or in the 
compact reddened absorption-line object Hawaii 167 (Cowie 
et al.[T994l . A counterpart to Al is SDSS 0127+01 14 (Hall et 
al. |2002 ), which is a strongly reddened mini LoBAL and has 
a Mg II configuration similar to the one observed in Al. In 
Hall et al. ( 2002 1 there is also a comparison between the dered- 
dened spectrum of SDSS 0127+0114 and a SDSS composite 
spectrum (their Fig. 12), which shows the strength of the ab- 
sorption. We performed a similar comparison by blueshifting 
the continuum-subtracted BL Lac spectrum by 25 A and scal- 
ing it to match the wings of the absorbed Mg II emission in Al . 
The comparison is shown in Fig. [5] As can be seen, the spectra 
match very well. The resulting FWHM would now be ~ 40 A 
instead of ~ 100 A. 

We investigated the morphology of A 1 on our I-band image 
in order to test for the presence of an AGN. The same analysis 
as for PKS 0537—441 was performed on the PSF-subtracted 
image (Fig. ^ right). According to our fits, an AGN+galaxy 
model is not strongly preferred over a pure galaxy model to 
describe the morphology of Al. Our best fit was an AGN+disk 
model with m core = 22.58, m^isk = 21.07 and r c = 0"61 (x 2 = 
1.17), whereas we got m disk = 20.9 and r c = 0"45 (x 2 = 1-36) 
for a pure disk model and nideVauc = 20.59 and r c = / .'64 (x 2 = 
1 .24) for a pure de Vaucouleurs model, respectively. We could 
not obtain a good fit with an AGN+de Vaucouleurs model. On 
the archival HST snapshot image of PKS 0537-441 Al and 
A2 are very faint, but Al appears more compact than A2. Thus 
there is no direct proof for a compact object in the center of Al, 
but this may be taken as further evidence for a (reddened) mini 
LoBAL QSO. 

If our favored explanation for the PKS 0537—441 - Al sys- 
tem is correct, this system would form a Quasar pair. Based 
on statistical arguments Kochanek et al. {T999 1 and Mortlock 
et al. ( 1999) have shown that the majority of wide separation 
(3"< A9 < 10") Quasar pairs are rather binary Quasars than 
gravitational lenses. Our system with a projected separation of 
3"9 and with a velocity difference of ~ 870 km/s would thus 
fall into the category of the rare class of binary Quasars. This 
is supported by the fact that the PKS 0537-441 - Al system is 
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a 2 R pair (both components detected in the optical, but only 
one in the radio; Kochanek et al. 1999), has the most extreme 
magnitude difference of all QSO pairs and differs strongly in its 
optical appearance. Given the relative small separation in pro- 
jected distance and relative velocity (the latter of which may 
even be smaller), this system may finally merge and form a su- 
permassive binary black hole (Begelman et al. 1980 1. 

We finally remark that the PKS 0537-441 - Al system 
would now already be the third binary QSO involving a BAL 
QSO (including LBQS 0103-275 3 (Jun kkarinen et al. l200!t 
and SDSS 0300+0048 (Hall et al. 120021 1. This is a high frac- 
tion of BAL QSOs amongst Quasar binaries with respect to the 
relative abundance of BAL QSOs. Curiously, all of these three 
Quasar binaries are at similar redshifts (z = 0.834, 0.858, z = 
0.892, 0.894, and z = 0.885, 0.893 for LBQS 0103-2753, SDSS 
0300+0048 and PKS 0537-441, respectively). 

The analysis of the cluster environment of PKS 0537—441 
indicates that the BL Lac is located in a cluster as rich as Abell 
type 0-1. The most comprehensive study on BL Lac cluster en- 
vironments up to z = 0.65 until now has been presented by 
Wurtz et al. dl997t . They found that most BL Lacs reside in 
a relatively poor cluster environment of Abell type < with an 
indication of evolution (richer environments) towards higher 
redshifts. Similar conclusions have been drawn by Fried et al. 
( 1993 1 for BL Lacs from the 1 Jy sample at z < 0.8. Thus our 
observations are in agreement with the observed trend. This is 
also supported by the detection of four galaxies at very similar 
redshifts (within A z < 0.002) as the BL Lac. However, al- 
though the objects for the MOS-setups were selected serendip- 
itously, an even higher number (10 galaxies, see Fig. |5J was 
found to be in the redshift range z = 0.9 - 1. According to 
the magnitude-redshift relation by Eales (1985 1 discussed in 
section 13. 1 .21 m„ (1) would only differ by 0.15mag between z 
= 0.893 and z = 0.95 (without k-correction). We can not rule 
out that the enhanced galaxy density is due to a cluster at even 
higher redshift (z - 0.95). Therefore PKS 0537-441 might 
either be projected onto a galaxy cluster at higher redshift or 
even be part of a large-scale structure with an extension to- 
wards PKS 0537-441. 

5. Summary 

We have presented a deep I-band image of PKS 0537—441 
as well as multi-object spectroscopy in the field of 
PKS 0537-441. Our main results are: 

1) We could not detect a host galaxy of PKS 0537-441 
or any galaxy along the line-of-sight to the BL Lac down to 
a detection limit of mi = 19.0. None of the companion galax- 
ies to PKS 0537—441 is a lensed background galaxy. Thus, we 
find neither evidence that the properties of PKS 0537—441 are 
affected by gravitational lensing nor that a potential group of 
galaxies surrounding PKS 0537—441 act as a lensing system 
itself. 

2) The nearby environment of PKS 0537—441 is very com- 
plex. The three major findings are: i) We detected extended 
off-nuclear [O II] emission in PKS 0537—441. This is most 
likely due to photoionization by the active nucleus, although 
an attractive alternative - jet-cloud interaction by the counter- 



jet in PKS 0537—441 - can not be ruled out. ii) One compan- 
ion to PKS 0537—441 has a very broad, asymmetric absorbed 
Mg II line at z = 0.885. This source is most likely a QSO, possi- 
bly a mini low-ionization BAL QSO. If this can be confirmed, 
PKS 0537—441 would be the first BL Lac being a member of 
a binary QSO. iii) We further found two companion galaxies 
within a few arcsec at similar redshifts as the BL Lac (z = 0.892 
and 0.895 vs. z = 0.893), which may trigger the strong activity 
in this source. One of the companion galaxies is at the location, 
where the radio jet of PKS 0537—441 turns over. 
3) PKS 0537—441 seems to be embedded in a galaxy clus- 
ter as rich as Abell type 0-1. This is supported by the detec- 
tion of four galaxies at similar redshifts as the BL Lac (Az < 
0.002). Since we found serendipitously an even higher num- 
ber of galaxies between z = 0.9 - 1, the BL Lac may also be 
projected onto a galaxy cluster at somewhat higher redshift (z 
~ 0.95) or be part of a large-scale structure with an extension 
towards PKS 0537-441. 

Clearly, the interpretation of the most interesting findings 
is difficult and somewhat speculative. High-resolution spec- 
troscopy of the extended [O II] emission line region (possibly 
also in [O III]) in PKS 0537-441 with the VLT are required to 
unambiguously discriminate between the photoionization and 
jet-cloud interaction model (e.g. by looking for kinematic sig- 
natures of a counterjet). A deep high-resolution VLA image of 
PKS 0537—441 would be useful to search for a counterjet di- 
rectly. High-resolution spectroscopy in the blue and moderate- 
resolution spectroscopy in the far-red part of the optical spec- 
trum or even in the NIR is necessary to confirm the presence of 
a (mini low-ionization BAL) QSO and to establish the presence 
of a binary QSO incorporating the BL Lac. 
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